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Abstract. We calculate the free energy, spatial string tension and Polyakov loop of the gluon plasma using the dilaton 
potential of Ref. [1] in the dilaton-gravity theory of AdS/QCD. The free energy is computed from the Black Hole solutions 
of the Einstein equations in two ways: first, from the Bekenstein-Hawking proportionality of the entropy with the area of the 
horizon, and secondly from the Page-Hawking computation of the free energy. The finite temperature behaviour of the spatial 
string tension and Polyakov loop follow from the corresponding string theory in AdS$ . Comparison with lattice data is made. 
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INTRODUCTION 

The duality of string theory with 5D gravity is nowadays 
a powerful tool to study the strong coupling properties 
of gauges theories. In conformal AdS$ the metric has a 
horizon in the bulk space at rj = %l 2 T where I is the 
size of the AdS-space, and entropy scale like s<* r\°zT 3 . 
To extend this duality to SU(N C ) Yang-Mills theory, one 
of the first tasks is to control the breaking of conformal 
invariance. In this work we will study the 5D dilaton- 
gravity model introduced in Ref. [2], with the action 



y = 



d 5 xy/^G [R- -d^tpd^tp - V (0) 




(1) 



87TG5 JdM 

The dilaton potential V (<p) is related to the j3 -function in 
a one-to-one relation. We assume the parameterization 
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which was proposed in Ref. [1] for the computation of 
the QQ potential at T = within this model. Eq. (2) 
has the standard behavior in the UV, j3(a) ~ —fiooc 2 — 
J3ia 3 + . . . for a -C a, and generates confinement in the 
IR. The optimum values to reproduce the QQ potential at 
T = are b 2 = 2.3, a = 0.45 and i = 4.389GeV- ! [1]. 

The model of Eq. (1) has two types of solutions. 
The thermal gas solution corresponds to the confined 
phase, and the black hole (BH) solution characterizes 
the deconfined phase and it has a horizon in the bulk 
coordinate at z = Zh [2]. The BH metric takes the form 
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where z — £ 2 /r and f(z) has the properties f(zh) — and 
f'izh) = —4-nT. Starting from the ultraviolet expansion 
of the j3 -function, one can solve analytically the equa- 
tions of motion in the UV [2, 3, 4, 5]. One arrives at the 
following expansion for the metric factor up to ff(a 2 ) 

{a) ^*\\-tfoa + ^(22ft-9h)cA . (4) 



FREE ENERGY 

The phase transition from the confined glueball gas to 
the deconfined gluon plasma can be understood as fol- 
lows. At low temperatures the physics is determined by 
the thermal gas metric with / = 1 . At a minimal temper- 
ature the BH solution appears with a horizon at Zh and 
/ 1, and becomes the preferred solution at the phase 
transition. So, the phase transition arises from the com- 
petition of the free energies computed in both metrics. 

We show two different ways to get the thermodynam- 
ics of the 5D dilaton-gravity model of Eq. (1). One pos- 
sibility is to use the Bekenstein-Hawking entropy for- 
mula which establishes the proportionality between the 
entropy and the area of the horizon of the BH. It reads: 
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where a/, = a(z/,). In the second equality we have used 
Eq. (4). From the entropy one can derive the pressure 
by solving dp(T)/dT = s(T). Another possibility is to 
compute the free energy from the Einstein-Hilbert action. 
One has to regularize the action by introducing a cutoff 
z = £ in the integral over the on shell action Eq. (1), c.f. 
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FIGURE 1. Free energy density as a function of temperature. 
We show as a full (red) line the numerical result obtained 
starting from the Bekenstein-Hawking entropy formula, c.f. 
Eq. (5). Blue points correspond to the result using Eq. (6). We 
include up to fi f (o&) in the r.h.s. of Eq. (7) to compute G. 
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FIGURE 2. T/^a~ s as a function of T (in units of T c ). We 
show as points the lattice data for SU(3) taken from Ref. [7]. 
The colored curves represent the analytical result from the 
holographic model, c.f. Eq. (10), and the black solid line refers 
to the full numerical result of Eq. (9). 



Ref. [2]. The free energy is computed as the difference 
between the free energy of the BH solution and that of 
the thermal gas solution. The result is [2, 5] 
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where Cf = 4nTb 3 (zi t ) corresponds to an enthalpy con- 
tribution in JT, and G= ({TtF^) t - (Tri^ v ) ) 
is the gluon condensate up to normalization factors. This 
method to compute & is more involved giving the fact 
that the computation of the gluon condensate G is rather 
subtle. The finite temperature solution of the Einstein 
equations differs from the zero temperature one at ff{z A )'. 
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In Eq. (7) care has to be taken to keep track of log- 
arithmic effects which are not usually considered, c.f. 
Ref. [2]. These higher order terms in (Xq affect appre- 
ciably the extraction of G from a comparison of the nu- 
merical solutions of b and bo. Only when at least the 
order ^(05o) is taken into account in Eq. (7), one gets 
good agreement of the thermodynamic quantities with 
the numerical results from the Bekenstein-Hawking en- 
tropy formula Eq. (5). We show in Fig. 1 the free en- 
ergy obtained by using the Bekenstein-Hawking entropy 
formula, Eq. (5), and the Einstein-Hilbert action from 
Eq. (6), using NNLO terms in Oo in Eq. (7) to compute 
G as a function of T. These plots correspond to a numer- 
ical solution of the equations of motion. See Ref. [5] for 
details. One recognizes in this figure the first order phase 
transition at T c = 273 MeV for zero flavours, which is 
quite close to lattice simulations. The upper branch in 
Fig. 1 represents the small black holes. 



SPATIAL STRING TENSION AND 
POLYAKOV LOOP 

The spatial string tension is very useful to test AdS/QCD 
models. It is non-vanishing even in the deconfined phase, 
and it gives useful information about the non perturbative 
features of high temperature QCD. The computation of 
the correlation function of rectangular Wilson loops in 
the (x,y) plane leads to a potential between quark and 
antiquarks which behaves linearly at large distances, i.e. 
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For details on the computation we refer to e.g. Ref. [6]. 
The spatial string tension takes the following form: 

a s (T) = ^at /3 b 2 (z h ). (9) 

From Eqs. (4) and (9), the ultraviolet expansion leads to 
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We show in Fig. 2 a plot of T/yfo~ s as a function of 
temperature including several orders in Eq. (10), and the 
full numerical computation from Eq. (9). One can see 
that the model reproduces very well the lattice data in the 
regime 1.10 < T/T c < 4.5. A good fit to the lattice data 
taken from Ref. [7] is obtained by using l s = 1 .94 GeV~ 1 . 

The vacuum expectation value of the Polyakov loop 
serves as an order parameter for the deconfinement tran- 
sition in gluodynamics. The correlation function of two 



Polyakov loops taken in the large distance limit leads 
to the vacuum expectation value of one single Polyakov 
loop squared. This means that the Polyakov loop is re- 
lated to the free energy of a single quark F q as 
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One can compute the Polyakov loop from the Nambu- 
Goto (NG) action of a string hanging down from a quark 
on the boundary into the bulk. The NG action reads [5, 8] 
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where we have reguralized it by substracting the action 
of the thermal gas solution up to a cutoff z c - The renor- 
malized Polyakov loop then writes 
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We show in Fig. 3 the behavior of Lr as a function of T 
computed numerically from Eqs. (12)-(13). We perform 
a fit to lattice data by considering the string length l s 
and z c as free parameters. The best fit in the regime T c < 
T < 10T C leads to l s = 2.36GeV~\ z c = 0.43 GeV -1 . 
Our approach fits Lr very well without a dimension two 
condensate, since a dimension two operator would show 
up in the UV expansion of the metric near z = 0, Eq. (7). 
This does not exclude that a good fit to the data exists 
of the form -21ogL R ~ a + b(T c /T) 2 with a = -0.23, 
b = 1.60, in accordance with Ref. [9]. See also Ref. [8]. 

The UV asymptotic behavior of the Polyakov loop can 
be computed from Eqs. (4), (12) and (13). It reads 
r £ 2 



L R (T) =exp 
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We plot in Fig. 3 also the UV asymptotics given by 
Eq. (14) and the pQCD result up to 0{a%) [10, 1 1]. 



CONCLUSIONS 

We study the thermodynamics of the 5D dilaton-gravity 
model by using a parameterization for the dilaton poten- 
tial which was quite successful to reproduce the QQ po- 
tential at T = [1]. We compute analytically in the UV 
and numerically the entropy, free energy, spatial string 
tension and Polyakov loop. We demonstrate the numer- 
ical agreement between computations of the free energy 
from the Bekenstein-Hawking entropy formula and from 
the Einstein-Hilbert action. Both approaches lead to the 
same result, but the later is much more sensitive to nu- 
merical errors, and an accurate computation is only pos- 
sible when one takes care of including logarithmic ef- 
fects in an UV expansion of the scale factor at finite T. 
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FIGURE 3. Renormalized Polyakov loop as a function of T 
(in units of T c ). Full (black) line corresponds to the numerical 
computation of Eqs. (12)-(13). We show as points lattice data 
for SU(3) taken from Ref. [12]. The colored curves represent 
the analytical result from the holographic model, c.f. Eq. (14). 
We also plot the result from pQCD up to 0{a?f) given in 
Refs. [10, 11]. 
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